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Isgur-Wise Function on the Lattice 
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We review our method and numerical results for calculation of the Isgur-Wise function on the lattice. We present 
a discussion of the systematic errors. Using recent experimental results, we find V c t = 0.044 ± 0.005 ± .007. 



1. Method 

A very active subfield in high energy physics 
recently is the study of hadrons with heavy-light 
quark content A major effort has been spent 
in calculating the Isgur-Wise function, which, 
once it is determined, can be widely used in cal- 
culations of heavy meson decay (b — > c) pro- 
cesses. After an initial exploration ||, calcula- 
tions of the Isgur-Wise function on the lattice 
H have quickly obtained interesting results which 
can be directly compared with experimental data 
and can be used to determine one of the elements 
of the CKM matrix, V c t, in the Standard Model. 

For calculating the Isgur-Wise function, £(v ■ 
v'), on the lattice, we have proposed Q to use 
the flavor symmetry of the heavy quark effective 
theory (HQEFT) g and measure the D -> D 
elastic scattering matrix element 

< D v i\cj„c\D v >= m D C cc (p)£(v-v'; fx)(v+v') v ,(1) 

where mo is D meson mass, v and v' are four- 
velocities of the initial and final D mesons. The 
constant C cc (p) represents the QCD renormaliza- 
tion effect from the heavy quark scale to a light 
scale p. The calculation was performed in the 
quenched approximation using Wilson fermions. 
Both light and heavy quarks are treated as prop- 
agating. For details of the numerical simulation, 
refer to Refs. 



From the lattice point of view, calculating the 
elastic scattering matrix element has significant 
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advantages. In comparison to the B — > D pro- 
cess, the elastic process on the lattice has much 
less noise and therefore has smaller statistical er- 
rors. Furthermore, because of the exactly known 
value 



< D v \cry„c\D v >= 2m_o, 



(2) 



at the "zero recoil" point v' = v, the lattice ar- 
tifacts that are independent of momentum can 
be removed without ambiguity using Eq. (0) as 
normalization condition for lattice data A 
similar strategy for B — > D decay would have 
introduced an extra (unknown) 0(1/itiq) correc- 
tion. Therefore, not surprisingly, the most accu- 
rate data obtained so far on the lattice are from 
D — > D elastic scattering However, inelas- 

tic processes, such as B — > D and B — > D* , can 
be valuable consistency checks fi|,|| . 

2. Systematic Errors 

For analysis of the systematic errors, let us con- 
sider the slope, p 2 , of the Isgur-Wise function at 
y = v ■ v' = 0. A fit of the lattice data p| to the 



relativistic harmonic oscillator model |12 



£(») = 



y + l 



exp 



y+l 



(3) 



gives p 2 NR = 1.41(19). For a model independent 
determination of p 2 , one may choose to fit £(y) 
near y = 1 



£(y) = l-p 2 (y-l), 



(4) 
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and obtain |§ p 2 = 1.24(26). All the fits have 
taken account of the correlations between data 
points using covariance matrices. There are sev- 
eral potential sources of systematic corrections: 
quenching, scaling violation, light quark mass 
m q dependence, finite volume effect, heavy quark 
mass tuq dependence. 

Quenching. The error due to quenching is the 
most difficult to quantify. Although the effect is 
expected to be small if a scale such as /„• is set 
to the physical value (we use Ref. || to set the 
scale with /„.), a systematic numerical study is 
still lacking. We will not give an assessment on 
the quenching effect here. 

Scaling violation. Since by using the normal- 
ization condition Eq. (^|) all the momentum in- 
dependent lattice artifacts are removed and the 
remaining scaling violations are proportional to 
y — 1. Therefore, we expect the residual scaling 
violations to be small. A fit to data at j3 = 6.3 
and P = 6.0 found a difference of 13% for p 2 NR . 

A direct check on the Euclidean invariance on 
the lattice is to measure the ratio of the form 
factors /_//_)_. This ratio was found small and 
consistent with zero within large errors . 

Light quark mass (m q ) dependence. Our lattice 
data for £(y) are presented with m q set to the 
strange quark mass, m q — m s . These data are 
directly relevant to processes such as B s — > D s , 
B s — > D*. For B — > D, they have to be extrap- 
olated to the "chiral limit" m q = m Ur d,. An in- 
spection shows that the linear size of the physical 
volume is in the range of (lOOMeV)" 1 [|, there- 
fore, at m q < m s the finite size effect becomes 
important and contaminates the m q dependence. 

To estimate m q dependence we therefore use 
data obtained on the largest physical volume 
(24 3 x 39 lattice at /? = 6.0). We first estimate the 
shift in p 2 NR from m q to m q ' with both m q , m q i in 
the range of m s . Then this shift in p 2 NR is extrap- 
olated to the chiral limit. Using this procedure, 
we find p 2 NR decreases by 12% from m q = m s to 
m q = m u ^ || ■ It is interesting to note that the 
sign of this shift is opposite to the chiral perturba- 
tion prediction flCf] a nd in agreement to the bag 
model calculation |1 lf| . It is important to confirm 
this trend in the future with improved statistics. 



Finite volume effect. To estimate the finite vol- 
ume effect, we compare our data on 16 3 x 39 and 
24 3 x 39 lattices at = 6.0, n q = .154 (m q = m s ). 
There is a shift of 15% in p 2 NR . We expect that 
the finite size effect would be smaller at a heavier 
m q . Indeed, the shift in p 2 NR is reduced to 9% at 
Kq = .152. 

Heavy quark mass (tuq) dependence. Re- 
cent lattice calculations indicated that the heavy 
quark symmetry begins to set in in the neigh- 
borhood of the charm mass. The leading 
1/itiq dependence agrees with the expectations 
of HQEFT. We refer to Ref. || for discussions of 
specific examples. Therefore, simulation results 
obtained at the charm mass range can be used 
and extrapolated to the heavy quark limit. For 
the Isgur-Wise function the leading order 1/mg 
correction should be ~ (y — 1)Aqcd/toq ]§]• It 
should be relatively small for current lattice cal- 
culations y — 1 < 0.2. Indeed, comparing p 2 NR at 
m Q ~ 1.6GeV and 2.3GeV, wc find 15% shift. 

Summary. Adding up the above items in 
quadrature, the total systematic correction be- 
comes 29%. We have 

p 2 NR = 1.41 ±.19 ±.41, (5) 

where the first error is statistical and the second 
is systematic error. For linear fit, we get 

p 2 = 1.24 ±.26 ±.36, (6) 

We should point out that this 29% system- 
atic error is probably an overestimate. Our fit in 
Eqs. (||) and (||) have been performed with data 
at all /3, lattice size, heavy quark mass values. 
Therefore, the combined systematic error is un- 
likely to be much larger than the statistical error 
(.19). Indeed, though we use it as an indication 
of the systematic errors, the shift in p 2 due to 
each item discussed above is not statistically sig- 
nificant. To get a better analysis of the system- 
atic errors, one needs more data points and bet- 
ter statistics. At this point, our discussion of the 
systematic errors should be taken primarily as a 
discussion on the methodology; the estimates ob- 
tained are only qualitative. 

A comparison with continuum model calcula- 
tions is given in Ref. |3j. Clearly, lattice result 
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has reached similar, if not better, numerical ac- 
curacy as the continuum models for p 2 . 

Our result for p 2 NR is also consistent with a re- 
cent lattice calculation by UKQCD Collaboration 
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Fig. 1. We obtain JTJ 



3. Extracting V c b 

Although p 2 is useful for comparison with con- 
tinuum model calculations, it is less useful for 
getting V c b from the experimental data. Around 
y = 1, experimental data have the lowest statis- 
tical precision Q. On the lattice, a model inde- 
pendent determination of p 2 also tends to have 
larger uncertainty because only a few data point 
close to y = 1 can be used. 



Figure 1. Comparison of present lattice data 
(crosses) with experimental data (open circles) 
[7]. The solid line is a fit to the lattice data ac- 
cording to Eq. (3). 



In ARGUS (and CLEO) experiments, what has 
been measured is with the most accurate 

data obtained in the range 1.1 < y < 1.5. In the 
lattice calculation, we have obtained in the 
range 1 < y < 1.2. Therefore, at least in the 
range l.l<y<1.2we can directly fit the exper- 
imental data with the lattice data with only one 
unknown parameter V c b- One such fit is shown in 



\V r 



cb 



TB 



1.53ps 



0.044 ± .005 ± .007, 



(7) 



where the first error is due to the statistical and 
systematic errors in the lattice calculation and 
the second error is from the experimental uncer- 
tainties. The errors on the lattice data essentially 
reflect the spread of £ over different j3, lattice size, 
and heavy quark mass values. 
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